Partial decoherence in mesoscopic systems 
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The coupling of a mesoscopic system with its environment usually causes total decoherence: at 
long times the reduced density matrix of the system evolves in time to a limit which is independent 
of its initial value, losing all the quantum information stored in its initial state. Under special 
circumstances, a subspace of the system's Hilbert space remains coherent, or "decoherence free", 
and the reduced density matrix approaches a non-trivial limit which contains information on its 
initial quantum state, despite the coupling to the environment. This situation is called "partial 
decoherence" . Here we find the conditions for partial decoherence for a mesoscopic system (with N 
quantum states) which is coupled to an environment. When the Hamiltonian of the system commutes 
with the total Hamiltonian, one has "adiabatic decoherence", which yields N — 1 time-independent 
combinations of the reduced density matrix elements. In the presence of a magnetic flux, one can 
measure circulating currents around loops in the system even at long times, and use them to retrieve 
information on the initial state. For N — 2, we demonstrate that partial decoherence can happen 
only under adiabatic decoherence conditions. However, for N > 2 we find partial decoherence even 
when the Hamiltonian of the system does not commute with the total Hamiltonian, and we obtain 
the general conditions for such non-adiabatic partial decoherence. For an electron moving on a ring, 
with N > 2 single-level quantum dots, non-adiabatic partial decoherence can arise only when the 
total flux through the ring vanishes (or equals an integer number of flux quanta), and therefore 
there is no asymptotic circulating current. 



PACS numbers: 03.67.-a, 05.40.-a, 03.65.Yz, 03.67.Lx 

I. INTRODUCTION 

Quantum computation operates on information stored 
in "qubits" , which are superpositions of two basic quan- 
tum states^ Clearly, quantum computation requires the 
stability of the quantum state stored on each qubit, and 
therefore it can be used only while this state remains 
coherent^ Interactions between qubits and their environ- 
ment, including input-output measurement devices, can 
cause decoherence which destroys the information stored 
in the qubits. Attempts to avoid decoherence have led to 
studies of decoherence-free subspaces, within which the 
quantum state remains protectedi^ The states in such 
subspaces are practically decoupled from the sources of 
the decoherence, due to symmetries of the system. The 
existence of such states means that the system has only 
partial decoherence. 

Consider a qubit which is based on two single-level 
quantum dots£ Denoting the states on each dot by |1) 
and 1 2), the initial state of the qubit is written as 

\ip ) = cosa|l) +e J7 sina|2) , (1) 

with two real parameters a and 7. In a recent paper ^ 
some of us considered the coupling of such a two- 
quantum-dot system to a general environment, and found 
conditions for having partial decoherence. An appropri- 
ate tuning of the dot energies and the coupling energies 
between the two dots yielded conditions under which one 



is able to retrieve the quantum information stored on the 
dot even after a long time, despite the decoherence. 

An alternative system involves coupling the two-dot 
qubit to a third dot, |3), which acts as a "control 
register" £ Only the latter dot is coupled to the envi- 
ronment. Such a ring of three dots was considered in 
proposals^ to reduce the decoherence of the single qubit 
state. Indeed, an analysis of the time evolution of the 
three-site ring, with a coupling of the 'control register' 
to an external fluctuator^ has found partial decoherence 
under appropriate conditions. 

The above two systems are special cases of a meso- 
scopic system, which contains N single-level quantum 
dots (or any other system which has N quantum states) . 
In this paper we generalize the discussion of Ref. 0, and 
study the conditions under which such a system can ex- 
hibit partial decoherence. To study decoherence, we con- 
sider the time evolution of the system in the presence 
of the environment. To concentrate on the state of the 
mesoscopic system, one traces the total density matrix 
of the system plus the environment over the states of the 
latter, ending up with the N x N reduced density matrix 
of the system itself, 

p(t) = Tr w [m)){*(t)\] , (2) 

where is the combined state of the system and the 
environment. 
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In this paper we consider the following Hamiltonian: 

% = H + H env + VU cnv , (3) 

where the first two terms describe the separate system 
and environment, and the last term describes the cou- 
pling between them. Here, the operator V acts on the 
system, while the operator U onv acts on the environ- 
ment. Examples of such couplings include the coupling 
of the system to a vibrational mode, e.g. VU cnv = 
£|l)(l|(b + &'), where w creates an excitation of this 
mode,— or to electrons on a neighboring single-electron 
transistor, VU cnv = ^|l)(l|c^c, where c 1 creates an elec- 
tron on the transistor— 

It is convenient to work with the eigenstates \i) of H , 
with eigen- energies if 

Ho\i) = ~e e \e, . (4) 

From now on, wc denote these eigenstates by £ or 
and the states localized on specific quantum dots by n or 
to. If the system is decoupled from the environment (i.e. 
V = 0), then the off-diagonal density matrix elements (in 
the basis of the \£)'s) exhibit Rabi oscillations, 

Pu>(t) = e i &>-*t)t put (p) (5) 

(we use h = 1). However, the coupling to the envi- 
ronment, reflected by the last term in Eq. ©, implies 
that transitions between environment states can reduce 
the coherence of the system states— Such decoherence 
usually causes the decay of all the off-diagonal reduced 
density matrix elements to zero. This decay is often at- 
tributed to random fluctuations of the phases of these 
matrix elements, and therefore it is called " dephasing 11 . 
Although this decay has been explicitl y d emonstrated in 
many special cases [see e.g. Refs. Ir^8l lldll2l - fl4l ]. we are 
not aware of its general proof for an arbitrary environ- 
ment. 

Without special tuning, the diagonal reduced density 
matrix elements also approach an asymptotic limit which 
is independent of the initial quantum state. For example, 
if the environment is a heat bath at temperature T, the 
system usually relaxes, or thermalizes, towards thermal 
equilibrium, p u ,(t -> oo) = S u ,e-^ kT y E/[er e "' /(fcT) ]- 
At T = this leaves only the ground state, while at 
T — > oo this yields the fully mixed state, p u {t — > oo) 
i Y , Again, although this behavior is found in many 
numerical or approximate calculations, we are not aware 
of any exact statements proving it. 

Below we discuss situations, in which some diagonal 
matrix elements do not evolve towards these trivial lim- 
its. Instead, they remain time-invariant, yielding par- 
tial decoherence. A special case of such partial decoher- 
ence has already been discussed in the literature— - — 
If [% ,V] = 0, then also [H a ,H] = 0, and the corre- 
sponding Hcisenberg operator H (t) becomes time in- 
dependent. This situation has been called " adiabatic 
decoherence" ; 12 i 13 or " dissipationless decoherence" — It 



has been show n 12 ' 13 that in this case all the diagonal re- 
duced density matrix elements p u remain time- invariant. 

In the present paper we discuss the more general case, 
of non-adiabatic partial decoherence, when ['Hq,V] ^ 0. 
Even without this commutator, we find that some parts 
of the initial state information can be retrieved from the 
reduced density matrix at any time. Adiabatic deco- 
herence then turns out to be a special case of our gen- 
eral discussion, and we show that this is the only case 
which yields partial decoherence for N = 2& However, 
for N > 2 one can also have other situations. 

The general conditions for partial decoherence are dis- 
cussed in Sec. II. Section III demonstrates these con- 
ditions for a few examples. Generally, we consider a 
tight-binding Hamiltonian, with "site" energies e n and 
"hopping" matrix elements J nm , 

^o = E e »H- E( J ™»H+ h - c -) • (6) 

n (nm) 

For the special case of one-dimensional rings, which con- 
tain N single-level quantum dots, this becomes 

^o = E[ £ «i")("i-( J -«+ii")( n + 1 i+ h - c -)] * ( 7 ) 

n 

where n = N + 1 is identical to n = 1. In some cases 
we shall assume that the ring is penetrated by a mag- 
netic flux $ (measured in units of 4> a /(2Tt) = fac/e). 
The Aharonov-Bohm effect adds phases to the hopping 
terms, J„.„ +1 = ^nin+i e ^ n ' n+1 > where is real and 

$ = Yln=i n+i- This flux generates a circulating cur- 
rent around the loop, which is discussed in Sec. IV. Ad- 
ditional comments on the results are given in Sec. V. 

II. PARTIAL DECOHERENCE 
A. General 

Generally, the time evolution of the reduced density 
matrix is given by 

p u ,(t) = Tr env (l\V(t))(*(t)\e') 

= TT cnv (e\e~ mt y(o))(*(oy m \e>) . (8) 

This expression becomes simple if the two energy eigen- 
states \£) and are both also eigenstates of V, 

V\£)=v t \t)- (9) 

In this special case, Eq. (JS) can be written asi2 

PiAt) = Tr onv [ e -^ t (£|vl/(0))(*(0)|ne^' t ] 

= Tr onv [e^' t e-^ t (£|vI/(0))(*(0)|O] , (10) 

where H e = e ( l + H cnv + v e U cnv acts only on the environ- 
ment states, and 1 is the unit operator there. For I = £' , 
one has e l ^'*e _ * w «* = 1, and thus 

Pu (t) = Tr env (£|f (0))<*(0)$ = p u (0) . (11) 
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Since 

p u = J2(£\n)p nm (m\l) , (12) 

nm 

Eq. (|11[) represents a combination of the reduced den- 
sity matrix elements (in the site representation) which 
remains constant at all times, preserving some of the ini- 
tial information. Furthermore, if v £ — v e , then Eq. (|10p 
reproduces Eq. ([5]), so that the whole Hilbert sub-space 
spanned by the states \£) and \£') remains entangled^ 
and decoherence free. 

The number of time-independent invariants thus de- 
pends on how many basis states \£) are also eigenstatcs 
of V. The maximal possible number is N. In that case, 
Eq. ([§]) holds for all £, which implies the commutation 
[H , V] = and thus adiabatic decoherence. It follows 
that Eq. (|TT1) also holds for all £. Since normalization 
requires the relation p u = 1 (wc allow only one elec- 
tron in the system, and this electron stays there forever) , 
we end up with N — 1 independent relations among the 
density matrix elements p nm {t), which remain invariant 
at all times. However, the number of time-invariants can 
have any value between and N — 1, as we demonstrate 
in the following sections. 

As explained above, adiabatic decoherence implies that 
all the diagonal reduced density matrix elements remain 
time- invariant, while all the off-diagonal elements usu- 
ally decay to zero. Therefore, this situation has been 
called "pure dephasing" . Even in this simpler case of 
decoherence, a decay of off-diagonal matrix elements 
has been proven only for specific models with specific 
approximations . 12 ' 13 It has also been found in our calcu- 
lations for non-adiabatic decoherence, where we replaced 
the environment by a fluctuator which is described by 
telegraph noised It would be useful to have a general 
analysis of this decay, including criteria for when it oc- 
curs. In this paper we proceed with the assumption that 
the off-diagonal elements (in the energy basis) indeed de- 
cay to zero. 

For non-adiabatic decoherence, the number of time- 
independent diagonal reduced density matrix elements 
may be smaller than N. When p u is not conserved, 
the system and the environment states remain entangled. 
This probably yields relaxation of these matrix elements 
into themalization. However, we are not aware of a gen- 
eral analysis of this process. Below we concentrate on 
the conserved diagonal elements. 



B. An alternative approach 

As we have shown, p u is time-invariant if the state \£) 
is an eigenstate of V, Eq. ©. One way to proceed is 
thus to find the states \£), and then apply Eqs. dU) to 
find conditions for partial decoherence. We employ this 
approach in Sec. IIIA below. Wc now present an alter- 
natice approach, which avoids the explicit identification 



of the states \£). This approach is shown to be more 
useful in the example presented in Sec. IIIB. 

Consider an hermitian operator 0, which acts only on 
the system, and assume that this operator commutes sep- 
arately with both Ji and V: 

[H 0) O] = 0, [V,O] = 0. (13) 

Since the Heisenberg equation of motion of this operator 
has the form 

d = i[H +VU env ,O} = , (14) 

it follows that O is independent of time, and therefore 
also in the Schrodingcr representation one has 

(tf(t)|0|*(t)) = (*(0)|O|*(0)) = const. , (15) 

independent of time. Taking the trace over the states of 
the environment then yields 

(*(t)\oMt)) =T*iMt))m)\o] 

= tr{Tr env [\nt))m)\0}} = tr[p(t)0] 

= P m n(t)O nm = Y Pmn{°)°nm > ( 16 ) 
ran ran 

where Tr and tr denote traces over all the states of the 
system+environment and over the states of the system 
only. The last equality gives a relation between the ele- 
ments of the reduced density matrix, p mn (t), which must 
hold at all times, implying partial decoherence. As we 
show below, the existence of an operator O which sat- 
isfies the above conditions depends on the relations be- 
tween the operators "H and V. 

We start with the condition [H Q , O] — 0. This implies 
that we can diagonalize "H and O simultaneously, i.e. 

0\£) = o t \£) , (17) 

with N arbitrary real eigenvalues {o £ }. The second re- 
quirement in Eq. (|13|) . [V,C] = 0, then implies 

(£'\V\£)(o i ,-o i ) = 0. (18) 

If all the off-diagonal elements of V are not zero, then all 
the o/s are equal to each other, and the only solution for 
O is proportional to the N x N unit matrix. Equation 
(|16p then implies that P(?(t) is a constant, but this 
sum is anyway equal to 1 due to the normalization. 

To have at least one diagonal element o f which differs 
from all the others, we need at least one column of off- 
diagonal matrix elements of V to vanish, 

(£'\V\£) = for a given I and for all £' ^ I . (19) 

This condition is equivalent to the equality (j9j), so that 
this particular state \£) is also an eigenstate of V. In 
this case, we can choose o t , = S u ,, i.e. O = \£}(£\, and 
Eq. (Til)|) implies that p u {t) = p u (0), as in Eq. (TIT]) . 
However, other eigenstates of W need not obey Eq. (|9|). 

So far, we have reproduced here the same results found 
in Sec. II A. However, below we show an example in which 
this approach is easier to use. 
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C. Adiabatic decoherence 

For TV = 2, V has only one off-diagonal matrix element, 
( 1 1 V 1 2) (in the basis of the system's energy eigenstates). 
Therefore, the states \t) become eigenstates of V only 
if (1|V|2) = 0. In this case, V becomes diagonal, and it 
commutes with T-L , which implies adiabatic decoherence. 
For TV = 2, one can write T-L = al+h er and V = cl+d-cr, 
where I is the 2x2 unit matrix, er is the vector of the 
three Pauli matrices and the coefficients are all real (for 
hcrmiticity) . The equation [H , V] = then implies that 
[b x d] = 0, and therefore - apart from a trivial constant 
shift in a (to the new value c/K), equivalent to a shift in 
the zero of the energy of the system - one must have V = 
KT-Lq, where K is & real number. For a given operator V, 
one can use gate voltages to tune the site and hopping 
energies of H (i.e. the three components of b, or the 
energy difference e 2 — e i an d the complex coupling J 12 in 
Eq. ©) in order to obey this condition^ 

For TV = 2, adiabatic decoherence requires the relation 
V = KH Q . Of course, a relation like V = KH Q implies 
adiabatic decoherence also for larger TV. However, for 
TV > 2 one can have adiabatic decoherence while V is not 
proportional to H (but still commutes with it). Again, 
we predict that all the diagonal reduced density matrix 
elements remain time- invariant. In both cases one might 
encounter situations in which some eigenvalues of V arc 
degenerate (see below). In that case, the corresponding 
off-diagonal matrix elements of p will exhibit Rabi oscil- 
lations, without any dephasing. The whole subspace of 
the degenerate states is then decoherence free. 



III. EXAMPLES OF NON- ADIABATIC 
PARTIAL DECOHERENCE 

The environment can couple to the system in many 
ways. Here we concentrate on the special cases in which 
the coupling of the system to the environment occurs 
only via a sub-system, of dimension M . For example, in 
the case TV = 3 we assume that only the 'control register' 
1 3) is coupled to the environment, hence M = 1. When 
the environment couples only to the sites n = TV — M + 
1, . . . , TV, this requirement becomes 

N 

(£'\V\£)= (£'\n)(n\V\m)(m\£) = , (20) 

n,m=N~M+l 

for a given £ and all I ^ £'. 



A. M = 1 

For M = 1, the environment couples only to the state 
| TV), and the right hand side of Eq. (|20[) contains only 
one term, with m — n = TV. This describes the exam- 
ples mentioned after Eq. In this case, V clearly 



does not commute with T~L . Since (iV| V]iV) 7^ 0, Eq. 
(|20l) holds only when we have (N\£) = 0. Indeed, if an 
eigenstate \£) of H Q is orthogonal to the state |TV), then 
V|^> = [JV><JV[V|JV><JV|^> = and the time evolution of 
this eigenstate is not affected by the environment. Such 
an eigenstate then belongs to the decoherence- free sub- 
space of the Hilbcrt space of the system. 

Consider now the Hamiltonian ©, and denote \£) = 
Sn^/vJ 7 " 1 )' with ip n = (n\£). The Schrodingcr equation 
for \£) becomes 

- e n)^n = -J2 J nm^m . ( 21 ) 
m 

for n = 1, . . . , TV. To achieve partial decoherence, we 
need states with ip N = 0. Therefore, the equation for 
n = N becomes 

E J NrJ>m = • (22) 
m 

The other N — 1 equations involve only the subsystem 
which contains the remaining N — 1 sites, without the 
bonds connecting them to the site N . Solving those 
equations yields the ^ TO 's, and substituting them into 
Eq. (f2"2")l then gives a relation among the J Nm 's which 
must be obeyed in order to achieve partial decoherence. 

Consider now the ring with N sites, in which only the 
state \N) couples to the environment. Using gauge in- 
variance, we place the whole flux on the bond (Nl), i.e. 
</Vn+i = Q$n.N- Equation (|2"2"|) becomes 

Jjge^i + ^Vi^v-i = • ( 23 ) 

The remaining JV — 1 equations, which now represent an 
open chain with N — 1 sites and with vanishing boundary 
conditions, have only real coefficients J„,„+i, and there- 
fore they always have solutions with real ip^s. Substi- 
tuting back into Eq. (|2"5]l . this implies that we can have 
partial decoherence only if z = e l * = ±1, i.e. if $ is 
equal to zero or to tt. Furthermore, given the VVi' s which 
result from these N — 1 equations, Eq. ([23]) also imposes 
a particular ratio between and J$ jy_ 1 , 

Jna/Jn^n-i - -^iV-iM ■ (24) 

For the specially symmetric case, with e n = and 

Jnli+i = J' t ne Hamiltonian of the open chain with N — 
1 sites is symmetric under n <-> TV — n. Therefore, its 
eigenstates should also obey ip n = ±ip N _ n with the two 
signs representing symmetric and anti-symmetric states. 
In particular, one has tp 1 — ± r tp N _ 1 . Having assumed 
that all the J's are equal, Eq. (fM]) implies that one has 
partial decoherence for the antisymmetric states if z = 1 
and for the symmetric states if z = — 1. 

Since adiabatic decoherence requires that [H Q , V] = 0, 
and since now (1|["H , V]|TV) = J t N ^ 0, it is clear 
that we cannot conserve all the diagonal /? w 's. There- 
fore, there is no partial decoherence for this case when 
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N = 2, where the only possible partial dccohcrence is adi- 
abatic. For N = 3, the symmetric and anti-symmetric 
cases relate to the 'bonding' and 'anti-bonding' states of 
the qubit, (|1) ± \2))/y/2. Partial decoherence is obtained 
for the state \tj = (|1) - z\2))/y/2, and therefore 

p u (t) = pn(t) + P22 W ~ 2zRep 12 (i) = const. (25) 



or, equivalently, 



P33(t) + 2zRcp 12 (t) = const 



(26) 



Indeed, these relations were found to hold in an explicit 
calculation of the time evolution of the N = 3 ring, using 
the telegraph noise model for the fluctuations in U env & 

The same result (f2l))) would follow from the method 
presented in Sec. IID. It is easy to check that for the 
symmetric ring with N = 3, the operator 



= |3)<3| + |1)(2| + |2)(1| 



(27) 



commutes with both H and V, and therefore Eq. p^|) 



yields Eq. |[25)>. 

For symmetric rings with 
within the open chain are ip n 
1,2, . . . , N — 1, with V n = 
there are (N - l)/2 (or N/2) 
[or (N — 2)/2] antisymmetric 
even). The diagonal elements 
symmetric (or antisymmetric) 
in time for z = —1 (or z = 1). 



larger TV, the solu tions 
= sm{-Kln/N)y/2/N, I = 

(-l)^- 1 Vi V -„- Therefore, 
symmetric and (N — l)/2 
states when N is odd (or 
of p corresponding to each 
state will remain constant 
These correspond to 



2 ^ • fTr£n\ 
N ^ S \~N~) 

n.rn 



7r£m\ 

11 



t (t)= const. (28) 



Writing O in the form 



O = 



°AA °AB 
°BA °BB 



the commutator with V is 



[V,0] 



-°AB V BB 
, BB°BA \ v BBi°Bb\ 



(33) 



(34) 



The solution of the equation v BB o BA = depends on the 
determinant oiv BB . In our case we assume Act[v BB ] ^ 0, 
and then o AB = 0, for any M. 

We next consider the commutator 



[n,o] 



AA> °AA 
O a A — O 



h 



BA U AA 



BB^BA 



AB U BB 



[h BB , o BB ] 



J AA^AB 



(35) 



where we set o AB =0. It is now convenient to diagonalizc 
separately the Hamiltonians of the two decoupled sub- 
spaces, 



L AA\ 



h 



BB I 



(36) 



The vanishing of the diagonal parts in Eq. (|35[) implies 
that these eigenstates can also be chosen to be eigenstates 



of O aa and (D BB , 



o A a = J^ ^^ ' 



J BB 



Equation (|34j) now requires the relation 

(P\VBB\P){Op-Op,) = 0, 



(38) 



For example, for N = 4 and 2 = 1 we have 
Pn + P33 - 2Rcp 13 = const. , 
while for z = — 1 we have two invariants, 



(29) 



P11 + 2 P22 + P33 + 2Rep 13 ± 2V2Rc(p 12 + p 23 ) = const. 

(30) 



B. M > 1 

For M > 1, it is easier to follow the method outlined 
in Sec. IIB. We now split the Hilbert space of the system 
into two subspaces, A and B, of dimensions N — M and 
M . In the basis of the sites, the bare initial Hamiltonian 
is 



H 



and the interaction is 



h 



AA 
l BA 



AB 
l BB 



V = 





v BB 



(31) 



(32) 



and the vanishing of the off-diagonal part of Eq. ([35 
requires the relation 



\^BA\ a ){°p -°a) = 



(39) 



The first of these relations is similar to Eq. p^|) . There- 
fore, it will be satisfied if there is at least one column of 
zero off-diagonal matrix elements in the matrix of v BB , 
or if all the o^'s are equal to each other. 

The second of these equations involves h BA . For the 
ring Hamiltonian ([T]). this matrix contains only two el- 
ements, namely J Nl and J N _ M+1 N _ M - Therefore, to 
have non-trivial values of o a or of Oa we need to satisfy 



'fi\N — M 



{0\N)J Nil {l\a) 
■l)J N -M + i <N -M(N-M\a) = 



. (40) 



Placing again the whole flux on the bond (iVl), the sep- 
arate wave functions for the open chains in the A and 
B subspaces are real, and therefore Eq. ([4T)]) can be sat- 



isfied only if z 



±1, together with a particular 



ratio between and J^_ M+1 n-m- 

For the specially symmetric case e n = and J, 



(0) 



J, one has ifi\N) 



n,ri+l 

±(8\N - M + I) and (all) 
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±(a\N — M). and therefore the final invariants will again 
involve symmetric and antisymmetric eigenstates within 
each subspace. Specifically, for N = 3 and M = 2 we 
find 

p n + 2zRcp 23 = const. (41) 
For N = 4 and M = 2 one has two invariants, 

Rc[p 12 ± p 34 ] = const. (42) 

IV. CIRCULATING CURRENTS 

As mentioned, in some cases the asymptotic state of 
the system maintains a circulating current, which can be 
used for retrieving information on its initial state. The 
site state \n) is an eigenstate of the position coordinate 
of the electron. Therefore, the position operator can be 
written as 

h = ^2n\n)(n\ . (43) 

a 

The current (or the velocity) operator is thus 

1 = dh/dt = i[H, h] = I s + I e , (44) 

where 

i s = i[n ,h] , i e = i[v,h]u em . (45) 

For the ring Hamiltonian ([7]), the first term becomes 

Is = -iJ2( J n,n+M(n + l\-h.c.) . (46) 
n 

The quantum average of this current is 

(/,) = Tr(pl s ) = 2 J2 Im(p„,„ + i J n+1 , n ) ■ (47) 

For simplicity, consider Eq. ([7]) in the fully symmetric 
case, e n = and J„°„ + i = J, and choose the gauge 
<j) n n+1 = = $/7V. The eigenstates of H are 

\Q =^2 e ^ en/N \ n )/^ > ( 48 ) 

n 

and i e = -2Jcos(0 + 2tt£/N), for £ = 0, 1, • ■ ■ , N — 1. 
Therefore, the conserved quantities are 

Pti(t) = PtM = E e 4(2 ^ /JV)(m -"V m J0)/A . (49) 

van 

Assuming that the off-diagonal matrix elements p u , de- 
cay to zero at long times, and the asymptotic reduced 
density matrix in the site representation becomes 

Pmn (t ^^)=J2 e l(27r£/Ar)( - m W(0)/iV . (50) 
e 



For N > 2, one may have degenerate states: Since 
~ H - e e , = 4Jsin[0 + tt(£ + t')/N] sin[7r(£ - t)/N], these 
two states are degenerate if cj> equals an integer multiple 
of 7r and ii £ + £' = N. In these cases, p u , (t) is also time- 
invariant. Specifically, for the adiabatic decoherence of 
the N = 3 symmetric ring with V = KH , the whole 
subspace formed by the states \£ = 1) and \£ = 2) (i.e. 
the states of the qubit) is decoherence free. However, if 
"i 7^ v 2 then only the diagonal matrix elements p xl and 
p 22 are conserved. 

For this special symmetric case, Eq. (|50l) yields 

(51) 

which is equivalent to the magnetic moment per site as- 
sociated with the flux through the ring. 

The second term in Eq. (|44[) . / e , represents current 
fluctuations due to the coupling to the environment. This 
term requires the average over the environment states of 
U CDV , and the commutator [V, n}. In the simplest adia- 
batic case, when V = KH , this will result in a simple 
shift of the prefactor J to Jn_ = J(l + K(U cnv )). This was 
the expression used in Ref . Jfj. In the case of coupling at 
a single site, Sec. IIIA, V is diagonal in the site basis, 
hence [V, n] = 0. In other cases, this term requires a 
separate analysis. 

For the special case ./V = 2, if we start with the initial 
state ([J), we have p u — [1 + sin(2a) cos(7r£ + 7)]/2, and 
thus p n = p 22 = 1/2, p 12 = p 21 = sin(2a)cos7/2. 
Therefore, (/) = J sin(2a) cos 7 sin <fiJ£- 

For all the non-adiabatic cases discussed above, we 
found partial decoherence only for real z = e l *. Since 
the asymptotic reduced density matrix is also real, such 
values of $ yield zero circulating current. However, this 
result is specific for single rings. As soon as the system 
contains more than one ring, the proof that e 1 * must be 
real breaks down, and one can still have some circulating 
currents. We expect to address this more general case in 
another publication. 



V. CONCLUSIONS 

We have shown that if a state \£) is simultaneously an 
eigenstate of both T-L and V then p u (t) is time- invariant, 
allowing the retrieval of some initial quantum informa- 
tion at any time, despite the decoherence. The number 
of such states \£) depends on the symmetry of the sys- 
tem, which can be modified by tuning the parameters 
in the system's Hamiltonian, and we have shown exam- 
ples in which one finds either adiabatic decoherence (all 
diagonal reduced matrix elements are time-invariant) or 
non-adiabatic decoherence. In the latter case, rings have 
an asymptotically vanishing circulating current. 

So far, we have concentrated on the ideal case, when 
the Hamiltonian can be tuned exactly so that one 
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achieves partial decoherence at all times. In Ref. \$ we 
also discuss small deviations from such exact tuning. For 
the special cases treated there we found that such devia- 
tions cause a very slow decay of the diagonal reduced den- 
sity matrix elements, which would remain time-invariant 
for the exact tuning. The rate of the decay depended 
on these deviations, and the amplitudes were equal to 
the initial values of these matrix elements. Thus, one 
could still extract information on the initial state from 
these amplitudes. We expect similar phenomena also for 
the more general case discussed here. This would result 
from an appropriate expansion of Eq. ([SJ in these devi- 
ations. Similar slow decays are expected to arise due to 
additional weak coupling terms, involving other environ- 
ments. We leave the analysis of such situations for future 
work. 

We hope that this paper will stimulate some experi- 



mental tests of our results. We also hope that it will 
generate some general discussion on what happens to the 
reduced matrix elements which are not time-invariant. 
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